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We investigate the relation between different three-dimensional reductions of the Bethe-Salpeter 
equation and the analytic structure of the resultant amplitudes in the energy plane. This correlation 
is studied for both the (j) 2 a interaction Lagrangian and the nN system with s-, u-, and t-channel 
pole diagrams as driving terms. We observe that the equal-time equation, which includes some 
£\J , of the three-body unitarity cuts, gives the best agreement with the Bethe-Salpeter result. This is 

followed by other 3-D approximations that have less of the analytic structure. 
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The interest in pion-nucleon (ttN) scattering has shifted in recent years from low energies to energies near and 
above the pion production threshold. This has been motivated by the growing importance of the low energy baryon 
^vi \ spectrum, and in particular the A* (1440) and A* (1535) in the study of QCD models. At these energies one needs 
to make use of a relativistic formulation of irN scattering, with the Bethe-Salpeter (BS) equation |l[ being a prime 
I candidate. The problem here is that the BS equation is four-dimensional, with a kernel that has a complex analytic 
r*"** ■ structure. This has inspired a number of three-dimensional reductions to the BS equation that have been employed 
in meson-exchange models of ttN and nucleon-nucleon (NN) scattering. In fact, it has been established that one can 
generate an infinity of such equations that satisfy two-body unitarity Q . The question that is often raised is which of 
these equations is the optimum one to be used for a given problem. To define such an optimum equation one needs 
, criteria for the definition of optimum. 

One such criterion proposed by Gross |!| stated that any relativistic two-body equation should have the correct 
'■^ \ one-body limit, i.e., in the limit as one of the masses goes to infinity the equation reduces to either the Dirac or 
Klein-Gordon equation. The Gross equation || achieves this limit by placing one of the particles on-mass-shell. This 
effectively reduces the BS equation to a 3-D equation, and in the ladder approximation has the correct one-body 
limit. The full BS equation does satisfy the one-body limit, but the ladder BS equation fails to satisfy this condition. 
This could be partly remedied in NN scattering by including the crossed as well as the ladder graphs in the kernel 
J> , of the BS equation ^, |5[ || . For irN scattering it is not clear as to how one may impose the one-body limit short of 
resorting to the inclusion of the crossed-ladder diagrams. Placing the pion on- shell § may be the optimal solution 
for iterating the u-channel pole diagrams, but is not as good for the t-channel pole diagrams [fi S|. 

An alternative scheme for developing relativistic 3-D equations that avoids the above ambiguities of the 3-D re- 
ductions, and is a natural extension of the non-relativistic two- body scattering problem, is based on a formulation 
of scattering theory using the generators of the Poincare group 0, [nj. This formulation, which introduces the 
interaction through one of the generators of the group, satisfies two-body unitarity and can be extended to include 
three-body unitarity. Recently, Fuda P4 [T^l has adapted this approach to the 7rA system within the framework of 
the coupled channel approach satisfying two-body unitarity only. To establish the connection with meson exchange 
models |L4|, the mass operator M in the interacting system is related to the mass operator in the free system (Mo) 
by the linear relation M = Mq + U in which U is the interaction. This interaction is defined on the Hilbert space 
of baryon (B) and meson-baryon (mB) states, allowing for the coupling between the two channels. To remove the 



energy dependence of the meson exchange potentials they resort to Okuba's |15 unitary transformation. The source 
of ambiguities in this approach are: (i) the choice for the Poincare generator in which the interaction is introduced, 
and (ii) the method by which the energy dependence is removed from the interaction. 

In the present analysis we propose to examine a third alternative to the reduction of the BS equation from a 
four-dimensional to a three-dimensional integral equation. Here we are motivated by trying to preserve as much as 
possible the unitarity content of the BS equation while preserving the two-body nature of the integral equation. This 
is partly a result of the observation that the low energy resonances observed in irN scattering are near the three-body 
threshold, and one may need more than just two-body unitarity. It also could form the starting point of a three-body 
formulation of ttN scattering without the complexity of the three-particle BS equation. 
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To examine the importance of the role of unitarity in the reduction of the BS equation, we examine a small class of 
three-dimensional equations that preserve some of the unitarity cuts in the BS equation. In particular, we consider 
the the Klein |L7| or "equal-time" (ET) equation || the Cohen (C) approximation jlSfl , the instantaneous (I) 
approximation 



>r equa 

|3, [2lJ used recently by Pascalutsa and Tjon for nN scattering [|2 lf| . and the Blankenbecler-Sugar 
(BbS) equation \22\. Since there are an infinite number of possible 3-D equations, there are many other equations 
that we could compare to the BS equation. However, we have chosen to consider only a small set of equations which 
differ in the analytic structure of the resultant amplitude in the energy or -^/s-plane. 

In Sec. II we first examine the analytic structure of the BS amplitude for <p4> scattering in the <p 2 a model by 
performing a pinch analysis p3[ on the relative energy integration in the kernel of the integral equation. We will find 
in this simple model that as we proceed from the BS equation to the ET, Cohen, and finally to the instantaneous 
equation and the BbS equati on, we are including less and less of the analytic structure present in the original BS 
amplitude. We then in Sec. Ill proceed to present the results of the same analysis for 7riV scattering, within the 
framework of s-, u- and t-channel exchange potentials. To illustrate the magnitude and effect of these thresholds for 
both cjxj) and ttN scattering, we present in Sec. |^ numerical results for both models over the region covering the 
lowest three-body thresholds. Finally, in Sec. we give some concluding remarks regarding the 3-D approximations 
to the BS equation. 



II. ANALYTIC STRUCTURE OF THE SCATTERING AMPLITUDE 



Unlike two-body equations in three dimensions, the BS equation generates unitarity thresholds by the pinching 
of relative energy integration contour by the poles and branch points of the integrand of the integral equation. To 
illustrate this, we consider a simple model in which the interaction Lagrangian is gcj> 2 a, where g is the strength of the 
interaction, while (f> and a are the two scalar fields in the model. We will further assume that the kernel of the BS 
equation for (jxp scattering is approximated by one cr-exchange. The motivation for this approximation is to simplify 
the analysis while maintaining the main features of NN scattering in terms of i-channel exchange, and nN scattering 
in terms of s-, u- and t-channel exchanges. 

The BS equation for (fxp scattering, after partial wave expansion, is a two-dimensional integral equation of the form 



T t {q,q';s) = V t {q,q' ]S )-i dq'J / dq"V e (q, q"; s) G(q" , s) T e (q" , q'; s) , 



(1) 



where q = (qo,q) is the zero component and magnitude of the space component of the relative momentum, while 
P = ki + &2 = {\/s, 0) is the total four-momentum in the center of mass. Although there are several definitions for the 
relative four-momentum pj| , the most convenient, from the point of view of the Wick rotation p5| , is q = \{k\ — k,2). 
For scalar particles, the two-body propagator in the center of mass is 



G(q,s) = 



go 



1—1 r 



qo 



-1 -1 



Ei 



(2) 



where E q = \J q 2 + to 2 , with m the mass of the </> particle. The partial wave potential due to the exchange of a a 
particle of mass \x is given by 



Ve(q,q';s) = 



0o - q'o) 2 -q 2 - q' 2 - 



(2tt) 3 " l V -2qq' 
Here Qi is the Legendre function of the second kind, and for I — this potential takes the form 



(3) 



V (q,q';s) 



9 



2(2tt) s 



log 



(go - gp)- - u 2 q+q , + 
(go-go) 2 - u\_ ql + 



(4) 



where ui v = yfp 2 + . This form exhibits the logarithmic branch points of the potential. We observe that the 
potential is independent of the total energy squared, s. 

The analytic structure of the partial wave amplitude in the energy {^/s) plane is determined by the pinching of 
the q'fy integration contour by the singularities of the integrand 26, ^7j. These singularities arise from the potential 
Vi(q,q"] s), the two-body propagator G(q",s), and the scattering amplitude Ti(q",q';s). The two-body propagator 
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has four poles in the ft -plane. These poles, labelled 1) . . .4), are located at 

1) ft = --yfs - E q „ + ie ; 2) ft = \V~s~ E q » + ie ; 

3) ft = -\^Ts + E q n-ie ; 4) ft = -y/s + E q „ - ie , (5) 

and are illustrated in Fig. |l|. As we increase y/s > the poles 1) and 3) move to the left in the g '-plane, while poles 
2) and 4) move to the right. 

To examine the conditions under which the ft integration contour gets pinched by the poles in the propagator, 
we consider Fig. [j]. Here we have to ensure that the two poles that pinch the contour of integration come from the 
Feynman propagators of two different particles. As a result poles 1) and 4) ( or 2) and 3) ) can pinch the contour, 
while 1) and 3) ( or 2) and 4) ) cannot pinch the contour as they belong to the same particle propagator. This 
results from the fact that while the Feynman propagator for a single particle can be decomposed into a positive and 
a negative energy component, these components cannot pinch the contour of integration as they cannot generate any 
new thresholds. We therefore have, for y/s > 0, the condition for the poles 2) and 3) to pinch the contour to be 



1 _ 1 

-y/s-E q „ +ie = -- 

or 



E q „ +ie = --y/s + E q/ , -ie (6) 



s/s = 2 E q „ - 2ie . (7) 

Here we note that if we deform the q" contour of integration in Eq. ([!]), e.g. q" — ► g"e - ^, then E q n will acquire a 
negative imaginary part and as a result the ft integration contour will not get pinched. This is true for all values 
of q" except q" — 0. Therefore we may conclude that the pinching of the ft contour generates a branch point at 
y/s = 2m, which is the threshold for two-body scattering. 

The same analysis can be carried through for y/s < 0. In this case the contour of integration in the ft variable is 
pinched by the poles 1) and 4), with the condition for the pinch being 

y/s = -2E q » + 2ie . (8) 

This will generate a branch point at y/s — —2m. This result establishes the fact that if we take both branch points 
resulting from the pinching of the integration contour by the poles of the two-body propagator, we get an amplitude 
that is analytic in s. In this way charge conjugation symmetry is preserved pgfl . However, if we only include the 
branch point at y/s = 2m, then we have only a right hand cut in the y^-plane and the resultant amplitude is a 
function of y/s, i.e., it satisfies two-body unitarity, but violates charge conjugation symmetry. 

The potential Ve(q, q"; s) gives rise to two pairs of logarithmic branch points in the g '-plane. These branch points, 
labelled 5) ... 8), are located at 



5) ft =qo- Uq+q" + ie ; 6) ft = g - ojq- q >> + ie 
j q+q n - ie ; 8) ft = g + 0J q - q " 



7) ft = go + Uq+q>> - ie ] 8) ft = g + ui q - q " - ie , (9) 



with branch points 5) and 6) ( 7) and 8) ) being connected by a logarithmic branch cut. These are illustrated in 
Fig. @. Since the potential is due to the exchange of a scalar particle, the potential, before partial wave expansion, is 



position of the propagator into positive 
and each will give rise to a separate 



proportional to the Feynman propagator for that scalar. Here again the deconr 
and negative energy components gives rise to the two logarithmic cuts in Fig. 
term in the multiple scattering series. 

We now can examine the pinching of the contour of integration in the g '-plane by a combination of a propagator 
pole and a potential branch point. Take for example the pole 2) with the branch point 8). The condition for the 
pinch is 

E q » +ie = g + uJ q - q >> - ie (10) 

or 

1 _ 

-Vs = <?o + E q » + u q - q " - 2te = p . (11) 

Here again contour deformation, e.g. q" — > q"e~ 1 ^, will stop us from pinching the contour since both E q n and ui q - q " 
acquire a finite imaginary part. The exception again is the point q" = 0. In this case the value of q" for which we get 
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the pinching of the contour depends on q. To determine the position of the branch point in the -y/s-plane, we need to 
find the minimum value of po or y/s with respect to q" for which we have a pinch, i.e., 

dpo _ jf_ _ (q - q") = 

dq" E q n LU q - q » 

This equation is used to determine q" in terms of q, and gives 

„ m 

q = — : — q ■ 

We now can establish that the off- mass- shell scattering amplitude has a branch cut in the qo-plane at p3f 

qo = -Vs - y/q 2 + (m + W + 2 *e ■ ( 12 ) 

This gives the singularity in the relative energy q' ' for the off-mass-shell amplitude that is being integrated over in 
Eq. §■ 

Following the same procedure, we examine the pinching of the contour by the singularities 3) and 6) for y/s > 0, 
while for y/s < the contour is pinched by either 1) and 8), or 4) and 6). As a result the off-mass-shell amplitude in 
the integral has four possible branch cuts in the g^'-plane. These are at 

9) q , ' = ^yfs-£ q n+2ie; 10) q' ' = ~yfs + £ q „ - 2ie (13) 

for y/s > 0, and for y/s < the branch cuts are given by 

11) q^ = ^ + £ q „-2ie; 12) <fi = -£,,, + 2ie , (14) 

where £ q = y/ q 2 + (m + fj,) 2 . In Fig. ^| we illustrate the positions of the corresponding branch points in the Qg-plane. 
At this stage we should recall that the singularities 9) and 10) are in different terms in the multiple scattering series, 
and therefore cannot pinch the contour of integration as may be inferred from Fig. |^. 

Since the singularities 9) . . . 12) are in the off-mass-shell BS amplitude, they are present in the integrand of the 
integral equation, Eq. (|l|). These singularities in the amplitude can, in conjunction with the propagator poles or 
potential branch points, pinch the q' ' integration contour and generate further singularities in the amplitude on the 
left hand side of Eq. (0). 

Let us first consider the pinching of the contour by the amplitude singularities and the propagator poles. For y/s > 
we now can have the integration contour pinched by the singularities 2) and 10) ( or 3) and 9) ). The condition for 
the pinching between the pole 2) and the singularity 10) is 



1 r- 1 

-Vs-E q n +ie= -- 

or 



-y/a-Egu +ie= --y/s + £ q » - 2ie (15) 



y/a = E q „ + £ q n - Me . (16) 

This pinching of the contour generates a singularity for q" = at \fs = 2m + /i. This is the three-body threshold 
for <r production. We can proceed in the same manner to examine the other possible pinching of the contour in the 
^Q-plane by the poles of the propagator and the singularities of the off-mass-shell amplitude. These will generate 
thresholds in the y^s-plane at 

y/a = ±(2m + y) . (17) 

This gives us one of the thresholds needed for the equation to satisfy three-body unitarity. At this stage our amplitude 
is a function of s and therefore satisfies the requirement of charge conjugation symmetry pi] ]. 

We can also pinch the q'^ integration contour by a singularity of the amplitude and a branch point of the potential. 
For example, if we pinch the contour of integration with singularities 6) and 10), we get a branch cut in the off-shell 
amplitude at qo — ^y/s + \/q 2 + (m + 2/i) 2 — 3ie. This in turn gives rise to a branch point in the off-mass-shell 
amplitude in the integral part of Eq. ([!]), and with the help of the propagator pole 2), will pinch the integration 
contour to generate the four-body threshold at y/s = 2m + 2[i. 
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One can continue to include higher singularities into the off-mass-shell amplitude by considering the pinching of the 
contour of integration by singularities of the amplitude and those of the potential. These off-mass-shell singularities 
when included in the integral part of Eq. ([!]) can, with the help of the propagator pole, pinch the contour in g ' , and 
generate further thresholds in the y^-plane. These thresholds are for two-, three-, ... a production and are located at 

s=(2m + n/x) 2 for n=l,2, .... (18) 

In this way we have demonstrated that the BS amplitude has not only the two-body threshold, but all ro-body 
thresholds. Unfortunately, in the ladder approximation with undressed propagators and vertices, the Bethe-Salpeter 
amplitude does not give the full contribution to n-body unitarity. 

We now turn to three-dimensional reductions of the BS equation (also known as quasi-potential equations), and 
illustrate how they have a smaller subset of the thresholds and/or off-mass-shell singularities than the BS amplitude. 
As a first approximation we consider the equal-time equation. This was originally suggested by Klein jlTj and has 
more recently been revisited for the (f> 2 a model [ jl8[ and the nucleon-nucleon interaction jfjj, where the problem has 
been restated in a more formal manner which allows one to consistently improve upon the 3-D results to ultimately 
achieve the BS result. In this case the 3-D equation is given by |l8| 

T 1 =K l +K 1 {G)T 1 , (19) 

where the 3-D amplitude T\ and potential K\ are defined in terms of integrals over the relative energy go, i.e., 

Ti = (G)- l {GTG) (Gy 1 

Ki = (G)~ 1 (GV G) (G)- 1 , (20) 
where the brackets in the above expressions are defined as 

(A) = dq dq' A(g ,q;g ,q';s) . (21) 



In this case the thresholds are generated by the double integral over the relative energies go and q' in the initial and 
final state in (GV G). Here, (see Appendix |a|) for ^/s > the q' integration will generate: (i) the elastic threshold 
at t/s — ±2m by the pinching of the two poles of the propagator to the right of V in Eq. (EOh; (ii) the pinching of 
the propagator poles with the potential branch points will generate the off-mass-shell singularities 9) and 10) in the 
go-plane for ^/s > 0. For y/s < 0, the singularities 11) and 12) are generated in the go-plane. These singularities, with 
the help of the poles from the propagator on the left of V in Eq. (|20|), will generate branch points at ^fs = ±(2m + /i) 
(see Appendix |A|) . 

Thus the ET equation includes, in addition to two-particle unitarity, the contribution to three-body unitarity from 
the cr-exchange diagram. Unfortunately, as was the case with the BS equation, these thresholds are only part of the 
three-body unitarity cuts that should be present. Additional contributions to three-body unitarity come from the 
dressing of the cj) propagators, the (j> 2 a vertex and from higher-order contributions to the potential such as the crossed 
box diagram. 

In the Cohen ||l9f approximation it is assumed that the amplitude in Eq. (jl]) does not depend on the relative energy 
g '. As a result, we can perform the g ' integration, i.e., 

dq'JVt(q,q";s)G(q",s) . (22) 

In this case the only pinching of the g ' integration contour is: (i) between two propagator poles, and (ii) between a 
propagator pole and a potential branch point. Here, in addition to the two-body threshold at y/s = ±2m generated 
by the pinching of the contour by the propagator poles, we have singularities in the off-mass-shell amplitude resulting 
from the pinching of the propagator poles and the potential branch points. These correspond to the singularities 9) 
and 10), 11) and 12). This introduces an inconsistency in the equation to the extent that the analytic properties 
of the amplitude as predicted by the integral part of the integral equation are not consistent with those assumed in 
the reduction of the dimensionality of the equation from four to three. To resolve this inconsistency it is assumed 
that the relative energy in the amplitude takes on its on-mass-shell value, i.e., go = 0. In that case the amplitude 
has, in addition to the two-body threshold, a branch cut at i/s = 2£ q ± 4ie and as a result a threshold in ^/s, at 
yfs = ±2(m + /i). This indicates that the Cohen amplitude has one of the four-particle thresholds, but no three-body 
thresholds. 

In the instantaneous approximation, the inconsistency has been removed by assuming a static potential and therefore 
fixing the value of the relative energy in the potential [^l| (e.g. go = <Zo =0). As a result the only singularity 
generated by the pinching of the g ' contour is by the poles of the two-body propagator. These generate the thresholds 
at y/s = ±2m. We should point out that for <jxf) scattering the Blankenbecler-Sugar equation is identical to the 
instantaneous approximation. 
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III. ANALYTIC STRUCTURE OF THE ttN AMPLITUDE 



We now turn to the analytic structure of the irN amplitude. Here we have three classes of diagrams that can 
contribute to the potential at the one particle exchange level: the s-, u-, and ^-channel pole diagrams. The s-channel 
pole diagrams do not contribute any singularities to the integral part of the BS equation, Eq. (|l|). As a result, the 
only thresholds generated by the s-channel pole diagrams in the potential are those resulting from the pinching of the 
integration contour by the irN propagator and amplitude singularities. These singularities are also present when the 
potential is due to u- and i-channel pole diagrams, and need not be considered until we examine the u- and f-channel 
pole diagrams. 

For the u- and t-channel pole diagrams, the potential takes the form 

v(q, i'\ P) = - ( J% q, ' lP l TT- > ( 23 ) 
[q + Qq) - m H + ie 

where £ = 1 for the it-channel pole diagrams, and —1 for the i-channel pole diagrams. Here ran is the mass 
of the exchanged particle. Thus for u-channel pole diagrams mn — JTi/v > "1a j ■ ■ ■ ; while for ^-channel diagrams 
fnu — * Mff = MpjMo-j""' The function F(q,q';P) depends on the spin of the exchanged hadron and the form of 
the coupling at the two vertices in the diagrams. Since the function F has no singularities that could take part in 
the pinching of the integration contour, we could assume F = 1 if we are willing to neglect form factors associated 
with the vertices, and we do not consider the values of the discontinuities across any cuts in the amplitude. This is 
equivalent to assuming all particles are scalars, and all physical thresholds are independent of the spin and isospin of 
the particles involved. Consequently, the BS equation maintains the form given in Eq. (|l|) with the contribution to 
the partial wave potential for u- and t-channel pole diagrams given for t — by 

V (q,q) = C, . 3 log<^ 7 — " 2 -, j-^ f-— \ . (24) 

2(2tt)^ L (g + C9o) - (? - C? ) - mjj + ie J 

The ttN propagator can be divided into two contributions, one from the positive energy component of the nucleon 
propagator, and the second from the negative energy component. We have in the center of mass 

G, w ( 90l q; S ) = G™ (9, S )A+(q) - G™ (g, S )A~(-q) , (25) 

where the energy projection operators are written in terms of Dirac spinors as 

A+(q) = y> r (qK(q) , (26a) 



A ~(q) = -I>r(qK(q) . (26b) 



Also 



GTn(Q,s) = ^ -—- -= ^— , (27a) 

G%{q,s) = ^ 1 r \ 2 7 —, (27b) 



with Eg = -\/q 2 + m 2 N and ui q = \J q 2, + /zf . In addition ajv and are functions of s such that ajv + a, = 1. The 
solution of the BS equation does not depend on the choice of ajv(s) and a w (s), so for convenience, in our discussion 
of thresholds generated by the BS equation, we use the simplest possibility: = a v = 1/2. 

We note that for G^n the number of poles is still four, and the pinching of the integration contour by the propagator 
poles generate the elastic thresholds at y/s = ±(totv + Mw)- Therefore it is only the pinching between the propagator 
poles and the potential branch points that needs to be considered. 



A. t-channel pole diagrams 

We first examine the thresholds generated by the i-channel pole diagrams. These are the result of pinching of the 
contour of integration in the q$ plane by the poles of the propagator and the branch points of the potential. Here 
we follow the same procedure adopted in the last section for <f><p scattering in the (f> 2 <j model, and state the results. 
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For y/s > 0, the pinching of the potential singularity with the propagator pole generates two branch points in the 
off-mass-shell amplitude at the relative energy 

Qo = T^Vs ± vY + (rn N + p H ) 2 T 2 i e . (28) 

These two branch points are also present in the off-mass-shell amplitude in the integral part of the integral equation 
Eq. ([!]), and with the help of either the propagator poles or the potential branch point, can pinch the integration 
contour to generate further branch points. In conjunction with the propagator poles the branch points in Eq. ( gg ) 
generate the threshold at 

yfs = (m N + p w + /is) . (29) 

These are thresholds for p and a production in a model in which the i-channel poles are represented by p and a 
exchanges. 

The pinching of the amplitude singularities in Eq. ( p8| ) and the potential branch points will, for yfs > 0, generate 
branch points in the off-mass-shell amplitude 

qo = T^Vs± vV + (m N + p H + p H ,) 2 =F 2ie . (30) 

Here, as in the case of scalar particles considered in Sec. II, the Bethe-Salpeter equation has all the thresholds for the 
production of the mesons present as i-channel pole diagrams in the potential. Thus the amplitude branch point in 
Eq. ( |30| ) with the help of the propagator pole will generate the four-body threshold at 

y/s = (m N + p^ + p H + m>) > (31) 

which is present in the off-mass-shell box diagrams given by G V* G Vt G. 



B. u-channel pole diagrams 

Turning to the u-channel pole diagrams, the branch points of the potential with the help of the propagator poles 
generate two off-mass-shell branch points at 

qo = T^Vs ± vV + ( m N + m H ) 2 T 2 i e . (32) 

These amplitude branch points in conjunction with the propagator pole generate three-body thresholds at 

\fs = 2rriN + run and ys = 2p v + ran ■ (33) 

To generate the four-body threshold we have to examine the three classes of off-mass-shell box diagrams corresponding 
to G V u G V u G,GV t GV u G and GV u GV t G. The first class gives rise to the threshold at 

\fs = (niN + P-k + tub + tub' ) • (34) 

This corresponds to first pinching the integration contour by the singularities in Eq. ([32]) and the branch point of the 
it-channel pole potential. This is to be followed by the pinching of the integration contour by the resultant branch 
point and the propagator pole. The second class of box diagrams GVtGV u G will generate thresholds at 

\G — (2m N + m H + Ph') and y/s — (2 p n + m H + Ph') ■ (35) 

Finally the third class of box diagram GV u GVtG gives four-body thresholds at 

\fs — (2m N + m H + Ph>) and ^fs = (2p n + m H + Pw) , (36) 

which are identical to those in GVtGV u G as expected. In the above run = tojv, wia, • • • and pu = p p , p a , ■ ■ ■ . We 
could continue this procedure to generate all the higher order thresholds in the Bethe-Salpeter equations as we did in 
Sec. II. 

The above analysis is for ^fs > 0, and similarly we could generate the mirror thresholds by examining the case 
y/s < 0. In this way we establish that the amplitude is a function of s rather than y/s. 

The point to be emphasised in establishing the positions of the thresholds in the s-plane is that these thresholds, 
included in the ladder BS equation, are not sufficient for the amplitude to satisfy three-, four- and higher-body 
unitarity. This means that we have some three-body thresholds, but we don't have three-body unitarity. Only two- 
body unitarity is included completely. To include three-body unitarity we would need to carry through a Faddeev 
decomposition of the amplitude, which would result in a set of coupled equations that are significantly more complex 
to solve g|. 
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C. Thresholds in the 3-D equations 

We now turn to the three-dimensional reduction of the BS equation for the ttN system. We consider the equal-time, 
Cohen, instantaneous, and the Blankenbecler-Sugar equations, each of which has different unitarity structure in the 
s-plane. The BbS equation is just one example of an infinite number of 3-D equations that are derived from the BS 
equation by replacing the irN propagator G^n with a new propagator containing a (5-function which fixes the relative 
energy. As a result, the only threshold included is the two-body unitarity cut. For more information about the 3-D 
equations, see Appendix 

The instantaneous and BbS equations generate thresholds at y/s — ±(rriN + /z^), however, the coupling to negative 
energy states is neglected in the BbS equation. The equal-time or Klein equation has, in addition to the two-body 
thresholds at y/s = ±{nij\r + ^n), the three-body thresholds as presented above for the ladder BS equation. This 
includes 

\/s = rriN + Hit + , Vs = 2niN + rriH , and \fs = 2/^ + tuh ■ (37) 

The Cohen approximation, in addition to the two-body thresholds at y/s — ±(m^ +/«„■), has thresholds resulting 
from the pinching of the contour of integration by the propagator poles and the potential branch points with q 
on-mass-shell, i.e., qo = (m 2 N — [i 2 ) / (2y/s) . For the i-channel pole diagram, we get for yfs > the branch point at 

y/s = m N + fa + 2/iff . (38) 

On the other hand, the it-channel pole diagram gives thresholds at 

yfs = — fj,^ + 3tojv + 2itih and yfs = — mjy + 3/^ + 2m# . (39) 

These are basically the four-body thresholds the Cohen approximation had in the scalar case considered in Sec. II. 
Here again the analysis is for yfs > 0, and can be extended to get the mirror thresholds for yfs < 0. 

IV. NUMERICAL RESULTS 

To examine the importance of the higher unitarity cuts as we proceed from the BS equation to the various 3-D 
reductions, we will present numerical results for the two models ((j) 2 a and nN) under consideration. Since thresholds 
do not depend on the spin and isospin of the fields involved, we have chosen the <fi 2 a model as an approximation 
to the nucleon-nucleon interaction with one pion exchange in the deuteron channel. There are no resonances in the 
system, and we expect the higher unitarity cuts to be of no great importance as is the case for the 3 Si NN channel. 
On the other hand, for the irN system the dominance of the A(1232) will allow us to examine the effect of the higher 
unitarity cuts on the width and position of a resonance. 

To emulate the NN system in the (f> 2 a model with a exchange as the potential, we have chosen the mass of the 
cj) (N) to be 1 GeV, while the a (n) has a mass of 0.15 GeV. For S-wave scattering we have adjusted the coupling 
constant so that the (f>(f> system has a bound state with a binding energy 2.2 MeV corresponding to the deuteron. In 
Fig. [| we present the phase shifts resulting from the solution of the BS equation as well as the equal-time, Cohen, 
and the instantaneous 3-D equations. Also included are the phase shifts resulting from the equal-time equation after 
excluding the negative energy component of the propagator (ET-NAP), i.e., we have excluded the left hand cut and 
hence violated charge conjugation symmetry. This situation corresponds to time-ordered perturbation theory p8[ , 
which has been used frequently in models of NN scattering (see, e.g., Ref. |3Cj) and nN scattering (see, e.g., Ref. J3 1[ ) . 
From the results in Fig. ^| we find, as expected, no great variation as one excludes higher unitarity cuts. On the other 
hand, the differences are consistent with the unitarity analysis of the last section in that the best approximation is 
the ET equation, followed by C, ET-NAP, and finally I. It is interesting to note that the inclusion of the three-body 
unitarity cut in the ET-NAP equation is more important than satisfying charge conjugation which is included in the 
I equation. 

A better test of the role of higher unitarity cuts is to examine the inelasticity for the different approximations to the 
BS equation, the results of which are shown in Fig. |^. In this case results for the I equation have not been included 
since the inelasticity is always unity, as this equations include the two-body unitarity cuts only. The agreement 
between the BS and ET equations is a result of the fact that both equations have the same three-body unitarity cuts. 
The Cohen equation does not have a three-body unitarity cut, and as a result, the inelasticity for this equation is 
less than one only for energies above the four-body unitarity threshold. At this stage the results for the BS and ET 
equations start to deviate as the latter does not include the four-body thresholds. Although we can observe differences 
between the equations in the inelasticity, in all cases the inelasticity is very small for this problem. 
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From the above analysis we may conclude that for energies close to the a production threshold, the ET equation 
handles the inelasticity much better than any of the other 3-D equations as it includes the same three-body unitarity 
cuts that the BS equation has. 

We now turn to the irN system, and consider a model in which the potential includes N and A s- and it-channel 
pole diagrams as well as a and p t-channcl pole diagrams. To solve the BS equation or one of the 3-D equations for the 
ttN system we need to introduce some form of regularization. This has been achieved by introducing two types 
of form factors for each of the vertices in the diagrams that constitute the potential. 

• Type I: Since the vertices in the potentials have three legs, we can write a form factor for each vertex as the 
product of three functions of the four-momenta of the three legs, i.e., 

fa/3-y (Pi ,Pp,Py) = fa (pI ) ffJ (Pp ) fy (P^) , (40) 

where a, f3, and 7 refer to the three legs of the vertex. The function f a (j>a) is taken to be 

'•<*-(^r- <4i) 

with n a = 1. Here m a is the mass of the particle, while A a is the cutoff mass that is a parameter of the 
potential. Since such form factors have poles, they can generate unphysical thresholds. To minimize the effect 
of these unphysical thresholds, we have constrained the cutoff masses A Q to ensure that these thresholds are at 
much higher energies than any of the physical thresholds we are examining. 

• Type II: To reduce the number of cutoff masses, and at the same time employ a form factor that is more 
commonly used, we have taken the form factor to depend on the four-momentum of the pion only, i.e., 

Mp 2 a ,p% Pl ) = U(pI) ■ (42) 

Here we present results for n = 4 only. 

The parameters of the BS potential and form factors have been adjusted to fit the ttN data up to a laboratory energy 
of 360 MeV |fl6|| . We use the BS parameters in the different 3-D equations, i.e., we do not refit the parameters using 
each of the scattering equations. This is so that the effects of the different 3-D reductions can be easily compared to 
the BS equation. 

In Fig. H we present the Sn and S31 phase shifts for the two types of form factors and the different 3-D approxi- 
mations to the BS equation. Also included are the results of the VPI SM95 j32| partial wave analysis. In this case 
it is not clear as to which 3-D equation gives the best approximation to the BS equation. This is partly due to the 
fact that the final amplitude is not dominated by any single one of the diagrams that contribute to the potential. In 
both the Sn and partial waves there are large cancellations between the various diagrams that contribute to the 
potential. 

To see a more pronounced difference between the different equations, we turn next to the P\\ partial wave. In this 
case the amplitude is dominated by the s-channel pole diagram, which is repulsive, as well as the w-channel A pole, 
which is attractive. In Fig. |t] we present our results for the two different types of form factor. Here it is clear that ET 
equation gives the best approximation to the BS equation, but the variation between the different equations is not 
significant provided we restrict our regularization to type II form factors that depend on the four-momentum of the 
pion only. For the type I form factor the BbS equation gives very poor results. At this stage we note that the BbS 
equation gives additional attraction when compared to the BS equation. 

We finally turn to the P33 partial wave which is dominated by just one term in the potential: the s-channel A pole 
diagram. Here (see Fig. ||) there is a dramatic change as one goes from the BS equation to the ET equation to the 
Cohen equation, and finally to the equations with just the two-body unitarity thresholds. The trend in the phase 
shifts suggests that the additional attraction generated by neglecting the higher unitarity cuts is slowly converting 
the resonance pole to a bound state pole. If it is the higher unitarity cuts that are the source of this difference, then 
these cuts would be even more important for higher-energy resonances in the irN system. This is a question that 
would require further investigation. 

For the case of type I form factors, tables of the numerical results are presented in Appendix 

We do not present results for the inelasticities in the irN system, because the inelasticity generated by the BS 
equation and the different 3-D equations is negligible for the energy range considered here. 
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V. CONCLUSION 



In the preceding analysis we have demonstrated that the Bcthc-Salpctcr equation in the ladder approximation has 
in addition to the two-body unitarity cut, a selected set of n-body unitarity cuts. In reducing the dimensionality 
of the BS equation from four to three it is possible to preserve some of the three- or four-body unitarity cuts and 
in this way preserve some of the features of the original equation. This could be important at energies close to the 
three-body threshold, as is the case for the irN system in the region of the lowest resonances. From the numerical 
results presented for the two models considered, we may conclude that the preservation of these higher unitarity 
cuts do in general result in an equation that is a better approximation to the BS equation. In particular, for the A 
resonance, the inclusion of these higher thresholds are of considerable importance. 

It is always possible to readjust the parameters of the potential within the framework of any one of these equations 
to fit the experimental data. This procedure will result in the coupling constants in the Lagrangian being different 
from those resulting from the use of the BS equation. As long as those coupling constants have little physical meaning 
one may justify using any scattering equation that preserves two-body unitarity. However, if we would like some 
consistency, e.g. between the irN and NN systems, then it may be important to reduce the approximations in the 
starting equation. More significant is the possibility of missing some physical mechanisms such as the contribution of 
three-body thresholds to resonances just above the pion production threshold. 

Although three-body unitarity is not completely included in the ladder approximation to the BS equation, one 
could add further contributions to three-body unitarity by dressing the nucleon in the ttN propagator. That will play 
an important role if the BS equation is to include coupling to other channels such as the it A channel. 
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APPENDIX A: THRESHOLDS IN THE EQUAL-TIME EQUATION 

To illustrate the thresholds generated in the equal-time formalism {| [l?], we consider the singularities of 

+00 

(GVG) = I d qo — * — * (Al) 

I (g-qo) -E 2 a +ie (^ + q ) -E?+ie 



+00 

dq' 



The singularity structure of (G V G) is determined by the pinching of the qo and q' contours by the poles in the 
integrand. For the q' integration, we have poles from both the propagators and the potential. As we have not partial 
wave expanded the potential we have at this stage only poles. However, upon partial wave expansion these poles will 
generate logarithmic branch points which wc will enumerate. 
In the <7q plane we have the four propagator poles: 

1) q' = ~Vs - E q > + ie ; 2) q = -y/s - E q , + ie ; 

3) q' = ~yfs + E q , -ie; 4) q' Q = ]pfs + E q , - ie . (A2) 

On the other hand the potential branch points in the ^g-plane are at 

5) q'o = qo - u q+q > + ie ; 6) q' Q = q - uj q - q , + ie ; 

7) q' Q = q + uj q+q , - ie; 8) q' Q = qa + u q - q > - ie . (A3) 
For y/s > 0, the pinching of the q' contour of integration is between the following singularities: 
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• Between two of the poles of the propagator, and in particular for yfs > 0, poles 2) and 3) pinch the contour 
provided 

= 2E q , - 2ie . (A4) 

In a similar manner for yfs < 0, poles 1) and 4) pinch the contour to generate a branch cut at yfs — —2E q + 2ie. 
These will generate thresholds at y/s = ±2m. 

• Between one of the poles of the propagator and a branch point of the potential. In particular singularities 3) 
and 6) pinch the contour if 

9) go = ~y/s + Eg> +u q - q , - 2ie . (A5) 

This branch point is going to be involved in the pinching of the go integration contour, and to that extent will 
not at this stage generate a threshold in the final amplitude. Although the go integration includes the on-shell 
point qo = 0, we do not expect this singularity to generate a threshold in the amplitude, since we can deform 
the go contour of integration to avoid the point go = 0. 

Similarly the propagator pole 2) and the branch point 8) pinch the q' contour provided 

10) qo = Ws-Eg> -uj q - ql +2ie . (A6) 

A similar set of branch points are generated for yfs < 0. In this case the propagator poles 1) and 4) (see Fig. |l|) 
pinch the q' integration contour to produce a branch point at yfs = —2m. On the other hand the pinching of the 
q' by the propagator pole 1) and the potential branch point 8) generates a branc h po int in the qo plane at the point 



qo = —\\fs — E q i — uj q - q i + 2ie. This is a reflection of the branch point in Eq. ( A6) in the -y^s-plane. In a similar 
manner, we get a reflection in the -^/s-plane of the branch point in Eq. ( |Aq ) when the pole 4) and the branch point 
6) pinch the q' contour. In this way the ET formulation maintains analyticity in s. 

We now turn to the qo integration. In this case the singularities of the integrand arise from the propagator poles 
and those generated by the q' integration. The propagator poles are at: 

11) qo = ~ E q + ie ; 12) q = iy^s - E q + ie ; 



2 

13) q = -^Vs + E q - ie ; 14) q = ^ + E q - ie . (A7) 



In this case the go contour can be pinched by the following singularities: 

• For yfs > 0, propagator poles 12) and 13) pinch the go contour to give the branch cut 

= 2E q - 2ie , (A8) 

while for < the poles 11) and 14) pinch the contour to generate a branch cut at yfs — —2E q + 2ie. These 
will generate the two-body threshold = ±2m. 

• The propagator pole 12) and the branch cut generated by the q' integration 9) can pinch the go contour. The 
condition for this pinch is 

y/a = E q + E' q + Lu q - q > - 3ie . (A9) 

Since this is in the kernel of the three-dimensional integral equation, the actual pinch takes place at q' = 0. As 
a result this will generate a branch point at yfs = 2m + /.t. 

In a similar manner the propagator pole 13) and the branch cut generated by the q' integration 10) can pinch 
the go contour. In this case the pinch takes place if 

= E q + Eg + ujq-q, ~ 3ie , (A10) 



which is identical to that in Eq. (A9), and will generate the three-body threshold at y/s = 2m + /x. A similar set 
of pinchings of the go contour takes place between the poles of the propagators and the branch points generated 
by the q' Q integration for y/s < to generate the three-body threshold at yfs = —(2m + fi). 

In this way we have established that the equal time equation generates an amplitude that has two- and three-body 
thresholds and is analytic in s. Unfortunately, not all possible three-body thresholds are included, and therefore the 
resultant equation does not satisfy three-body unitarity. 
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APPENDIX B: 3-D EQUATIONS FOR ttN SCATTERING 

For completeness, in this appendix we give explicit forms for the different 3-D equations discussed in the present 
work. Since some 3-D equations violate Lorentz invariance, for example the BbS equation, the solutions of such 
equations depend on the choice of the relative four-momentum. Here we use the commonly-used choice 

,s s + m 2 N — fi s + ix — m% 

a N{s) = , a n (s) = . (Bl) 

2s 2s 

In this section we make use of ttN — + itN amplitudes sandwiched between Dirac spinors. We use the following notation 
for the amplitudes: 

T au (q\q;s) = u(q')T(q\ q; s)u(q), (B2a) 
T^(q',q;s) = v(-q')T(q' , q; s)u(q), (B2b) 

and similarly for the potentials. With this notation, the BS equation has the form of two coupled equations for T uu 
and T° u 

T* u (q\q ] P)=V™ u (q\q;P)--^- [ ]T J d 4 q" > W ' (q' , q"; P)Gf N w " (q" ; P)T*" u (q" ,q; P), (B3) 



w"— u,v 



with w — u, v. The nN propagators G™ N W are defined in Eq. (|27^ 

Firstly we consider the Cohen equation, which has the form (before partial wave expansion) 

T fite (q , ,q; fl ) = V^(0,q , ;0,q; fl )-^u E [ dq" K^ w " (q' , q" ; s)^ (q" , q; s) , (B4) 

w"—u,v 

with the kernel given by 

/oo 
dq V*' w {0, q'; 90 , q; s)G™(q , q; s). (B5) 
- OO 

The relative energy integration can be carried out by making use of the Wick rotation [p5| . 
The remaining 3-D equations for 7riV scattering have the form 

(q,q;«)=A (q , q; s) - — — I 2^ dq K (q , q ; s)g 7TN (q;s)T (q,q;s). 

w"—u.v 



(B6) 



The two-body ttN propagators are denoted as g^ and g% v N , which contain the positive and negative energy compo- 
nents of the nucleon propagator, respectively. The 3-D potentials are denoted as K w w . 

To our knowledge the ET equation for ttN scattering has not been previously discussed in the literature. The two- 
body ttN propagator appearing in the ET equation is obtained from the propagator in the BS equation by integrating 
out the relative energy, i.e., 



/oo 
dgoG7rJv(<?o,q; s) 
-oo 

m N ( A+(q) A-(-q) 



2m- 



2E q w q \«/s — E q — L0 q + it s/s + E q + Lu q — ie / 
ee .g^(9;s)A + (q)+ 5 r w (g; S )A-(- q ). (B7) 

For the calculation of the ET potential, we will also need the inverse of g^N, which is given by 

5 ^(q; S )- 1 = ^^^([ A /i-^-^]A+(-q) + [^+^+ W9 ]A-(q) ] . (B8) 
2m E q \ J 

The potential for the ET equation is then 

/OO I* oo 

dq'o / dq Q G 7rN (q' Q ,q';s)V(q' Q ,q';qo,q;s)G^ N (q Q ,q;s) 
-oo J — oo 

xg^Niq-^y 1 . (B9) 
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Multiplying Eq. (B9) from the left and right by Dirac spinors, and making use of the orthonormality of Dirac spinors, 
we find 



K^(q',q; S ) = - (t> \ 2 ? " q q (V~s - t'E* - ^) (^- £E q - & q ) 
{2tt) z tun Tnjsf 

x dq'J dq a G^'(q' ,q'; S )V a ' w (q' ,q';qo,q;s)G^(qo,q; S ). (BIO) 



Here £' = +1 (—1) for w' = u (v), and £ = +1 (—1) for w = u (v). Again, the relative energy integrations above can 
be performed numerically by making use of the Wick rotation p5| , which can be applied in the same way as it is used 
in solving the BS equation jl6| . 

For the instantaneous equation, the are obtained by integrating out the relative energy of the ttN propagator 
appearing in the BS equation, and so we have 

E q Uq[y/S - L q - LU q + lfi\ 
91n{1\s) = ^ r r- - „ —y (Bllb) 

It is clear that the ET and I equations make use of the same irN propagator. 

The BbS equation is obtained by replacing G^n in the BS equation with the following two-body propagator 

G% s (q;P) = 2m / ~[a N f' + / 1 + m N ]SW[(a N P' + q) 2 ~ m 2 N }6W[( a7T P' - q) 2 - vl], (B12) 

/„ s' — s — it 

J s th 

where s t h = {w^n + ^) 2 , and P' = (Vs* /yfs)P. This results in the BS equation being reduced to a 3-D equation of 



the form of Eq. (B6), with the irN propagators given by 



hj q LO q \S - (hq + bJq) Z + lt\ 

9: v N (Q^) = 0. (B13b) 
Notice that the coupling to negative energy states is neglected, but has poles at both = E q + uj q and 

\fs = -E q - LO q . 

In the I equation the potential is assumed to be static, while for the BbS equation the relative energy is fixed at its 
on-shell value. Therefore the potentials for these two equations are given by 

K*'™(q', q; s) = ^'"(Q, q'; 0, q; a). (B14) 
The differences between the I and BbS equations lie in the form used for the 7riV propagator. 

APPENDIX C: TABULATED RESULTS 



Tables pi to IV show the numerical values of the nN phase shifts when type I form factors are used. 
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BS 


ET 


C 


I 


BbS 


50 


6.78 


6.30 


6.34 


7.34 


9.33 


100 


9.17 


7.84 


7.96 


9.18 


12.9 


150 


10.2 


8.05 


8.12 


8.76 


13.0 


200 


10.3 


7.42 


7.35 


7.39 


12.3 
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9.76 


6.21 


5.92 


5.40 


11.0 


300 


8.72 


4.62 


4.06 


2.95 


9.13 


350 


7.18 


2.74 


1.90 


0.191 


6.85 



TABLE I: The Sn phase shifts using the type I form factors as calculated by the BS equation and the different 3-D quasipotential 
equations. 
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-5.52 


-5.21 


100 


-8.73 


-5.96 


-9.24 


-9.34 


-8.45 


150 


-11.9 


-8.02 


-12.3 


-13.0 


-11.4 


200 


-15.0 


-9.96 


-15.1 


-16.5 


-14.4 


250 


-18.2 


-11.8 


-18.0 


-20.2 
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-13.6 
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-24.3 


-20.5 


350 
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Slab 


BS 


ET 


C 


I 


BbS 


50 


-1.09 


-1.21 


-1.50 


-1.07 


0.266 


100 


-1.53 


-1.40 


-2.46 


-2.34 


1.65 


150 


-0.50 


0.159 


-1.36 


-1.52 


4.60 


200 


2.45 


3.793 


2.03 


1.62 


10.9 


250 


7.75 


9.78 


8.26 


7.95 


23.6 


300 


15.7 


18.1 


17.7 


18.4 


45.2 


350 


25.9 


28.10 


20.7 


32.7 


70.9 



TABLE III: The Pn phase shifts using the type I form factors as calculated by the BS equation and the different 3-D 
quasipotential equations. 
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BS 


ET 


C 


I 


BbS 


50 


4.85 


4.61 


8.23 


174.0 


180.0 


100 


18.7 


16.7 


45.4 


171.0 


181.0 


150 


51.6 


42.2 


111.0 


169.0 


182.0 


200 


96.0 


80.5 


133.0 


167.0 


183.0 


250 


120.0 


108.0 


141.0 


165.0 


183.0 


300 


131.0 


121.0 


144.0 


164.0 


183.0 


350 


138.0 


131.0 


147.0 


163.0 


183.0 



TABLE IV: The P33 phase shifts using the type I form factors as calculated by the BS equation and the different 3-D 
quasipotential equations. 
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FIG. 1: The poles of the two-body propagator in the go'-pl anc - 
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FIG. 2: The branch points of the potential in the go'-pl anc - 
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FIG. 3: The branch points of the amplitude resulting from the pinching of the integration contour by the propagator poles and 
potential branch points in the g^-plane. The arrows indicate the direction in which the branch points move with increasing s. 
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FIG. 4: The S-wave phase shifts for 4><j) scattering with one cr-exchange as the potential using the BS, ET, ET-NAP, C, and I 
equations. 






FIG. 8: The P33 phase shifts for the type II (a) and type I (b) form factors. The data points are from 



